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Abstract 

Equations of motion for the microscopic number density /5(x, t) and the momentum density 
g(x, t) of a fluid have been obtained in the past from the corresponding Langevin equations repre- 
senting the dynamics of the fluid particles. In the present work we average these exact equations 
of microscopic dynamics over the local equilibrium distribution to obtain stochastic partial differ- 
ential equations for the coarse grained densities with smooth spatial and temporal dependence. In 
particular, we consider Dean's exact balance equation for the microscopic density of a system of in- 
teracting Brownian particles to obtain the basic equation of the dynamic density functional theory. 
In the thermally averaged equation for the coarse grained density p{x,t), the related dependence 
on the bare interaction potential in Dean's equation is converted to that on the corresponding 
direct correlation functions of the density functional theory. 
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The equations of fluctuating nonlinear hydrodynamics (FNH) [l|-|3| provide a very useful 
description of the time dependent behavior of various types of liquids, ranging from simple 
liquids J4I, solutionsjs, 6|, supercooled metastable liquids j? - 10 1, mixtures H], and complex 
liquids [12, isl. These equations primarily represent the time evolution of the respective 
coarse grained densities for the system. Experimentally measurable correlation functions of 
these densities are obtained from the equations of FNH by averaging the corresponding quan- 
tities with respect to the noise. The simplest example of a collective density is the number 
density p(r, t) of particles in a fluid. The latter is obtained as the average of the correspond- 
ing microscopic quantity, p(x, t)=p(x, t), where the hat over p indicates its dependence on 
the phase space coordinates of the fluid particles and the angular brackets denote an averag- 
ing over the appropriate ensemble. Dynamic equations for the microscopic densities (with a 
hat) are obtained exactly from the corresponding equations of motion of the individual fluid 
particles. For example, if the fluid particles follow reversible Newtonian dynamics (ND) the 
corresponding exact balance equations of hydrodynamics are known as Euler's equations 
which are also time reversible. For fluid particles with dissipative dynamics, the correspond- 
ing equations for the microscopic densities are also time irreversible. Apart from 
being exact representations of the corresponding microscopic dynamics, these exact equa- 
tions are formal and needs to be averaged over an appropriate ensemble for practical use. In 
this paper we obtain coarse grained forms of these exact equations by averaging them over a 
suitable local equilibrium distribution. The final equations are stochastic partial differential 
equations having smooth spatial and temporal dependence. The coefficients involved in the 
FNH equations are related to the thermodynamic properties of the fluid. 

We begin by considering a particle system in which the microscopic dynamics of the 
constituent particles is described by Smoluchowski equation [igI involving only the particle 
coordinates. Let the a-th particle (for a = 1, ...N) be of mass m and its position coordinates 
are given by {xq,}. A colloidal system with heavy particles in a solution is a typical example 
for such a system. In the over damped limit, the time dependence of the momentum of the 
particle is ignored and the stochastic equation of motion for the coordinate Xq, of the a-th 
particle is obtained as. 



7o- 



dt 



N 



5^V„f/( 
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The quantity 70 which has dimension of [MT~^] is chosen to be unity in the following. 
[/(xq, — x,^) is the interaction potential between particles a and v. The symbol with 
the Greek subscript a denotes the derivative operator with respect to the corresponding 
component of Xq. The stochastic part or the noise ^Q,(t) in the equation of motion ([T]) is 
white and gaussian. The noise correlation is given by, 

< >= 2kBT5^p5,,Sit - t') . (2) 

The collective density p(x, t) is defined in terms of phase space variables {'Xa}=TM(t) as: 

N 

/3(x,t) = ^5(x-x,(t)) . (3) 

a=l 

Since the number of particles is constant p(x, t) is a microscopically conserved quantity. 
Dean[l^, obtained the following equation of motion for p : 



-p(x, t) = A;bTV'p(x, t) + V ■ {p(x, t)V j rfx'f/(x - x )p(x', t)} + C(x, t) , 



(4) 



dt 

where V denotes the derivative operator with respect to components of x. The correlation 
of the random force denoted by C(^7^) is obtained as, 



C(x, t) = - J] V, ■ {(5(x - ^a{t))Ut)} = V ■ J] 5(x - Mt))Ut) . (5) 

a a 

C(x, t) is dependent on the phase space coordinates through Xq. Eqn. (jl]) is an exact repre- 
sentation of the stochastic dynamics described by eqns. ([I])-(l2l). Unlike the time reversible 
Euler equations, the Dean equation is dissipative since the corresponding microscopic dy- 
namics described by eqn. ([1]) is dissipative. We average eqn. (j4]) to obtain a partial 
differential equation for the corresponding coarse grained density p(x, t) which has smooth 
spatial and temporal dependence. This stochastic equation describes the time dependence 
of density fluctuations in a non-equilibrium state. The averaging should therefore be done 
over the corresponding non-equilibrium ensemble. We approximate the latter in terms of 
the local equilibrium ensemble. This is motivated from the fact that the hydrodynamic de- 
scription of the non-equilibrium fluid corresponds to the time regime in which it has reached 
a state of local equilibrium. At this stage the local densities are sufficient to describe the 
state of the system. The average of the microscopic densities (dependent on phase space 
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coordinates) over the local equilibrium distribution function /i.e obtains the corresponding 
hydrodynamic field with smooth spatial and temporal dependence. 

We note that the temperature T enters the exact eqn. (jlj) only through the first term on 
the RHS. This is a consequence of using the chain rule of the stochastic differential equations 
of the Ito calculus in obtaining eqn. Let a set of stochastic variables (i = 1, ..,m) 
satisfy the stochastic equation Xi = hi + gij^j with the correlation of the white noise being 
given by < ^i{t)^m{t') >= SimS{t — t'). Ito chain rulejl8| obtains the stochastic differential 
equation for the variable y{{xi}) in the form 

Edy . 1 d'^y 

i i,j,k 



With eqn. ([T]) for dynamics of x^, the matrix gij = ylkBTdij is diagonal. The first term on 
the RHS of eqn. (jlj) follows from the gij^s in the second term on the RHS of eqn. ([6]). This 
holds even for a system of noninteracting particles. We define next the averaging procedure 
with respect to the local equilibrium distribution /i.e. The latter is considered to have 
attained a stage when the system is at temperature T. This thermalization process is closely 
linked to the momentum distribution becoming the Maxwellian. So far the momentum 
current density has not been included in the formulation. The momentum current density 
g is the fiux for number density p(x, t) in the corresponding balance equation representing 
number conservation. 

^p(x, t) = -V. P"'^(^ - = -V.g(x, t) (7) 

a 

The probability function for the local equilibrium state is obtained in analogy with that 



of the equilibrium state [19|. We consider here the distribution function /i.e in terms of 
the corresponding local thermodynamic variables {/3(x, t), /i(x, t), t'(x, t)}. The distribution 
/i.e(r^;t) in the local rest frame, moving with the local velocity v(x, t), is obtained as 

/i.e(r'^; t) = Q-^exp ( - j drf3{r, t) [e'(r) - /x(r, t)p'(r)] ) . (8) 

where Q is the partition function Q = f dTfie(Tj^,t). The corresponding microscopically 
conserved densities of the energy and number of particles in the local frame are respec- 
tively denoted by e' and p'. The fluid has attained a state with a fixed temperature T 
{(3{r,t)=/3 = l/{kBT)) and the corresponding local equilibrium distribution is given by 



/i.e(r^;t) = Q ^exp^ — (5 H' — J drfi{r,t)p'{r) The equal time correlation of the mo- 
mentum density g'(x, t) in the local rest frame is now obtained as, 

(^i(xi,t)^j(x2,)f:))j^ = (5(xi - X2)5^jp(xi,i(;)A;Br (9) 

where p(x, t)=< /5(x, t) is the local density. 

Having defined the distribution function we now average the microscopic eqn. (jl]) over 
the local equilibrium ensemble to obtain an equation for the coarse grained particle density 
in the following form: 

ap(x, t) 



dt 



ksTV^pi^, t) + V-[j d^{VU{^ - x')}< p(x, t)p(x', t) >ij + C(x, t) . (10) 



The thermally-averaged equation (|TOl) has stochastic and regular parts. The stochastic part 
or the noise t) is obtained by coarse graining of the microscopic quantity t) defined 
in eqn. ([5]). Correlation of ( is understood as a combination of two steps. First, in ( we 
average the noise ^ of the microscopic equations of motion over the different configurations 
of the Brownian particles so as to obtain the coarse grained noise (. This is indicated with a 
subscript l.e on the angular brackets. Second, we correlate the coarse grained noise ( at two 
different space time points while the equilibrium temperature of the bath is maintained at T. 
We denote the latter with a subscript B on the angular brackets to indicate the averaging 
over the bath variables. The correlation of the noise ( in the coarse grained equation is 
obtained by interchanging the order of the two operations stated above. 



{C(x,«)C(x',(')> = ((C(x,«))_^(c(x',«')), 
- ((c(x.t)C(x',0), 

= 2A;BTVp(x,t)V'(5(x-x)5(t-t') . (11) 

V' denotes the derivative operator with respect to the components of x'. In the Markovian 
approximation of large separation between the characteristic time scales of the solute and 
that of the solvent variables, the noise correlation is independent of the coarse graining 
process. 

The first and second terms in the RHS of eqn. ( fTOil are in terms of slowly varying 
quantities. These are conveniently expressed in terms of the Liouville operator C For the 
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particle system £ is defined as 



N 



OL = l 1- 



(12) 



where the Hamiltonian H = K + V is a sum of kinetic and potential parts. The momentum 
dependent or the kinetic energy part is K = Yla=iPa/i'^^) interaction part V is 

a function of the particle coordinates only. The time rate of change of p(x, t) is expressed 
in terms of the operator C to obtain iCp = — V.g = p(x, t). Note however that if C acts 
on the momentum current density g, it does not obtain the corresponding time derivative. 
Using the definitions ( |T2|) and respectively of the Liouville operator C and the current 
g' in the local rest frame, we obtain 



N 

m^) = - 5^ V, Y,vVj{^ - x„) + y dx'{v.f/(x - x')}p(x, t)/)(x', t) 



(13) 



Next taking an average over the local equilibrium ensemble in the co-moving frame, we 
obtain 



;-z£^Kx,t))i,e = kBTVip{^,t) + j rfx'{Vif/(x-x')}<p(x,t)p(x',t) > 



l.e 



(14) 



Since £ involves derivative operators, integrating by parts and using the property i£H = 0, 
we obtain 



:-^a(x)) 



l.e 



rfr^,'(x)^i£exp -/3{H- / rfxV(x')p(x')} 



= P j dxV(x')(^:(x)z£p(x'))i., 

= -/?/' rfxV(x') (^:(x)V;.^;.(x')),^ = p(x)V.p(x) . (15) 

i 

In reaching the last equality we have used the key relation (|9]) for momentum correlation in 
the local equilibrium distribution which has thermalized at temperature T. 

The quantity on the RHS of eqn. ( fT5|) is a thermodynamic property. We link this to the 
Helmholtz free energy F expressed as a functional of the inhomogeneous density p(x). Using 
the equilibrium relation F — G = = —PV, where Q is the thermodynamic potential, we 



have in the density functional formalism 
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2l| 



F[p(x)] ^ fi[p(x)] + J rfxp(x)/x(x) . (16) 

Q[p] is a functional of the density obtained from the equivalent result of grand canonical 
ensemble partition function i7[p(r)] = — fc^TlnS. The density functional theory identifies 
the equilibrium density by minimizing the grand potential [SQ/Sp{'x)] = 0. Using the above 
relations it then follows that the corresponding Helmholtz free energy functional satisfies 
[5F[p]/5p(x)] = /i(x). Eqn. ( !TO|l then reduces with the help of eqns. (fT^ and (fT5|) to the 
form 

Pin Tap 

+ C ■ (17) 



pV— 

6p 



We have now put the bare diffusion constant Dq = 70/ (/3m) in the RHS above to keep the 



dimensional factor explicit 22 



23| . The free energy F[p] of the density functional theory is 



expressed ( in units of /c_bT) as a sum of two parts F[p] = Fid + -^ex, respectively denoting 
the non- interacting or ideal gas part -Fid[p] and the interaction part -Fex[p]- The ideal gas 
part is obtained as 



Fidlp] = ksT j dxp(x)[ln(A3p(x)) - 1] . (18) 

where Aq = h/ \/2/Kmk^ is the thermal wave length. The interaction part or the so called 
excess part -Fex[p] is generally expressed in a functional Taylor expansion in terms of direct 
correlation functions c^*)(xi, ...,Xj) at a density po- 



cW(xi, ...,Xi;po) 



5p(xi)...5p(xi) 



(19) 



P=PO 



The two point function for i = 2 is the Ornstein-Zernike direct correlation function c*^^-*. 

Going beyond the over damped limit considered above, the equations of motion of the 
fluid particles are obtained in terms of both the respective momentum and position coordi- 
nates. This is termed as the Fokker-Planck dynamics (FPD). For the particle system the 
momentum of the a-th particle is taken as XQ,(t) ( mass taken as unity) and its time rate 
of change is 



dt 



N 



Vv„t/(x,(t) -x^(t)) -7op„ + ^«(t), 



/3=1 



(20) 
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where 70 is a dissipative coefficient. Note that this model does not conserve momentum 
microscopically unless the dissipative coefficient 70 has a operator associated with it. 
Balance equations for the collective densities {p{'x,t), g{'x,t)} have been obtained by Naka- 
mura and Yoshimori[l5|. By averaging these microscopic balance equations, we obtain the 
corresponding equations of motion for the coarse grained variables {p(x, t), g(x, t)}. 



^p(x,t) + V-g(x,t) =0 
d 



(21) 



^^,(x, t) + 7o^7^(x, t) + / {J2 ^^(x - x„(t))} 

\ a 

+ j dx'{ViU{^ - x')}< p(x, t)p(x', t) >i ^ = 0i(x, t) . 



l.e 



(22) 



The noise ^j(x, t) in the generalized Langevin equation is obtained by coarse graining of the 
noise 9i[^,t) which is defined in terms of the noise ^q, in the micro-dynamic equations. 

N 

^,(x,t) = J]5(x-x,(t))e(t) . (23) 

Using the same approximations as discussed with eqn. (fTTi) . the correlation of the multi- 
plicative noise 6'j(x, t) is obtained as, 



(x, t)^,(x', t')) = 2fcBT7op(x, t)5(x - x')5(t - t') . 



(24) 



We evaluate the third term on the LHS of eqn. (!22|) by making a change of variables 
P« = Pa + ^v(r'Q,) and = r'^ in the locally moving frame with velocity v(r'Q,). Using the 
symmetry of the distribution function ([8]) for the locally moving frame under the p'— t- — p' 
transformation this term reduces to 

"5'i(x,t)fi'j(x,t)" 



A;BrVip(x,t) + V,- 



p(x,t) 



(25) 



Substituting the above result in eqn. ( 1221) . and using the relations (1141) -(1151). which also hold 
for the Fokker-Planck dynamics, we obtain the following stochastic equation for the coarse 
grained momentum density 



P 



+ pVi— +7ofi'i = fc'i 
6p 



(26) 



Eqn. (p6|) is very similar in form to the corresponding equations of fluctuating nonlinear 
hydrodynamics for the particles which follow reversible NDjs], [2^. There is however one 



important difference. Unlike in the case of the FPD, for the ND the dissipative term involves 
a 1/p nonlinearity jsj. On the other hand, in the FPD the noise correlation given by eqn. 
f l24p involves the density and hence the noise is multiplicative. The presence of the 1/p 
nonlinearity in the ND case is important with respect to the ergodicity-nonergodicity (ENE) 
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26f of dense liquids. The latter 
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28| . The ENE transition being 



transition of the self-consistent mode coupling theory 
has been extensively used for studying glassy relaxations 
referred to is a consequence of feedback effects from slowly decaying density fluctuations 
in a dense liquid. It has been established ^| that the 1/p nonlinearity in the generalized 
Langevin equation is essential in smoothing off this transition. For Newtonian dynamics 
of particles, the single equation fll7p for p has also been obtained 2J| by eliminating the 
current g from a field theoretic formulation of the problem. This however requires in the 
ND case applying the so called adiabatic approximation in which momentum fluctuations 
decay much faster than the density fluctuations. 

In summary, the present work fills up an important gap between a) a nonlinear diffusion 
equation for the coarse grained density field in terms of thermodynamic direct correlation 
functions and b) the exact representation of the Brownian dynamics of a system of particles 
(e.g. Dean's equation). Thus the single equation f lT7|) for p(x, t) is the coarse grained form of 
Dean's exact equation for the microscopic density p(x, t). The first equation has been used 
in a large number of works on fluids and is well known to the community as an ingredient 
of the dynamic density functional theory. The driving free energy for this coarse grained 
equation is obtained here in the standard form of density functional theory. The role of 
the bare interaction potential pU{r) in Dean's equation is replaced by the corresponding 
direct correlation function c(r). The averaged equation which forms the basis for field 
theoretic models for the dynamics of fluids, has often been considered (without proof) to be 
exact for interacting Brownian particles. The effects of nonlinearities in the coarse grained 
equation are studied with various analytic techniques which are available for calculating 
renormalized time correlation functions in terms of appropriate mode coupling contributions 
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These models often form the basis for studying dynamics heterogeneities and 



growing dynamic correlation lengths in supercooled liquids. 
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